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Pool  boiling  heat- transfer  measurements  were  made  using 
seven  Gewa-T  copper  tubes  in  R-113.  The  first  tube  (19-mm 
OD  *0.75  fins/mm  *0.25  nun  gap)  was  tested  under  three  condi¬ 
tions:  (a)  plain,  (b)  with  four  shrouds  and  (c)  with  up  to 

5  wires  (0.13  mm  diameter)  wrapped  in  the  inter-fin  trough 
area.  Verifying  prior  data,  the  shroud  with  longitudinal 
openings  of  60°  at  the  top  and  8.5°  at  the  bottom  gave  the 
best  performance  among  the  four  shrouds  tested.  This  shroud 
increased  the  boiling  heat-transfer  coefficient  by  253  per¬ 
cent  (over  the  smooth- tube  case)  at  a  heat  flux  of  10,000 
2 

W/m  while  it  showed  only  a  18  percent  increase  at  100,000 
2 

w/m  .  When  wire  wraps  were  provided,  in  all  cases,  the 
heat-transfer  coefficient  was  improved  at  all  heat  fluxes. 
The  use  of  three  wires  gave  the  best  performance  with  341- 

percent  increase  in  heat-transfer  coefficient  at  10,000 

2  2 
W/m  and  a  130  percent  increase  at  100,000  W/m  . 

The  next  three  Gewa-T  tubes  had  dimensions  of  25  mm 

OD  x0.75  fins/mm  with  fin-tip  gaps  of  0.15,  0.25  and  0.35  mm 

respectively.  The  last  three  Gewa-T  tubes  had  the  same 

dimensions  except  that  the  fin  density  was  1.02  fins/mm. 

For  each  of  these  fin  densities,  the  0.25  mm  fin-tip  gap 

produced  the  best  boiling  performance  at  all  heat  fluxes. 

Also,  for  a  given  fin-tip  gap,  the  boiling  performance 
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increased  with  increase  in  fin  density.  For  a  fin-tip  gap 
of  0.25  mm,  the  lower  fin-density  tube  produced  an  80-percent 
increase,  while  the  higher  fin  density  tube  produced  a  103- 
percent  increase  in  the  boiling  heat-transf er  coefficient. 
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Ra 

Pr 
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average  wall  superheat  (K) 

temperature  at  the  base  of  the  fins  (°C) 

fluid  saturation  temperature  (°C) 

average  wall  temperature  (°C) 

heat- transfer  rate  from  the  cartridge  heater 
(W) 

heat- trans f er  rate  through  unenhanced  end  (W) 
total  heat-loss  rate  (W) 

2 

heat  flux  from  the  Gewa-T  surface  (W/m  ) 

2 

average  heat- transfer  coefficient  (W/m  .K) 

thermal  conductivity  of  the  copper  (W/m.K) 

thermal  conductivity  of  Freon-113  liquid 
(W/m.K) 

volumetric  thermal  expansion  coefficient 
( 1/K) 

2 

kinematic  viscosity  (m  / s) 

2 

thermal  diffusivity  (m  / s) 
density  of  Freon-113  (kg/m^) 

Nussel  number 
Rayleigh  number 
Prandtl  number 

voltage  across  the  cartridge  heater  element 
(volts) 

voltage  drop  across  the  precision  resistor 
(volts) 

current  through  the  heater  element  (amps) 
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0.  D. 
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o 
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i 
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LU 


A_ 
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resistance  of  the  precision  resistor  (ohms) 

outside  diameter  of  the  Gewa-T  tube  (m) 

diameter  at  the  position  of  the  thermocouples 
(m) 

diameter  at  the  base  of  the  fins  (m) 

outside  diameter  of  the  unenhanced  ends  (m) 

inside  diameter  of  the  unenhanced  ends  (m) 

outside  wall  perimeter  of  the  unenhanced 
ends  (m) 

fin  height  (m) 

length  of  enhanced  test  section  (m) 
length  of  unenhanced  ends  (m) 
corrected  length  of  unenhanced  ends  (m) 
tube  thickness  of  unenhanced  ends  (m) 

2 

surface  area  at  the  base  of  the  fins  (m  ) 

cross-sectional  area  of  tube  wall  at  the 
unenhanced  ends  (m2) 

outside  surface  area  of  the  unenhanced  ends 
(m2 ) 

2 

gravitational  acceleration  (m/s  ) 
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excluded) ,  and  the  two  thermocouples  for  vapor  and  liquid 
temperature  measurements  were  calibrated  in  an  insulated, 
stainless-steel  bath.  The  tubes  were  completely  immersed 
into  the  bath,  which  was  provided  with  a  motor-driven  mixer. 

The  bath  temperature  was  measured  using  a  platinum-resistance 
thermometer  accurate  to  iO.OI  K. 

Since  the  temperature  measurements  during  this  investiga¬ 
tion  ranged  from  45°C  to  60 °C,  thermocouple  calibration  was 
performed  in  the  range  of  45°C  to  70°C.  To  begin  the  cali¬ 
bration  process,  hot  water  (at  about  75°C)  was  added  to  the 
insulated  container.  Following  proper  mixing,  all  thermocouple 
readings  were  automatically  recorded  through  the  data- 
acquisition  system.  The  bath  temperature  was  measured  using 
the  platinum- resistance  thermometer.  To  obtain  other  lower 
temperatures,  small  quantities  of  water  were  added  to  the 
bath,  allowing  5  minutes  for  temperature  stabilization  of  the 
test  section. 

The  discrepancy  (i.e.,  the  platinum- resistance  thermometer 
reading  minus  the  thermocouple  reading)  was  plotted  against 
the  thermocouple  reading  as  shown  in  Figure  4 .  The  second 
order  polynomial  curve  shown  in  this  figure  was  gem  raced 
using  the  data  from  all  seventeen  thermocouple j .  However, 
for  clarity,  only  three  thermocouples  were  selected  for  data 
points  from  each  tube.  It  is  evident  that  the  thermocouples 
in  both  tubes  show  similar  discrepancies,  the  thermocouple 
readings  have  only  about  ‘0.1  K  discrepancy  around  the 
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Plexiglas  lid,  which  had  been  fitted  with  a  rubber  O-rinq. 
Several  holes  were  drilled  in  the  lid  to  fix  in  place  the 
test  section  extension,  vent  tube,  connections  for  primary 
and  secondary  condensers,  and  two  more  thermocouples  to  sense 
both  the  bulk  liquid  and  vapor  temperatures.  The  Pyrex 
glass  vessel  with  the  attached  fittings  was  placed  on  the 
plate  heater,  which  served  to  maintain  the  liquid  at  satura¬ 
tion  temperature.  A  secondary  condenser  was  mounted  external 
to  the  vessel. 

C.  INSTRUMENTATION 

The  power  supply  was  connected  to  the  cartridge  heater 
through  a  voltage  regulator  and  then  to  a  Variac  scaled  from 
0  to  120  VAC.  This  arrangement  allowed  fine  adjustment  for 
a  stable  feed  voltage  during  the  runs.  A  precision  resistor 
of  2.031  ohms  was  connected  in  series  with  the  cartridge 
heater  to  provide  a  measure  of  the  current  flow  by  reading 
the  voltage  drop  across  this  element.  A  digital  voltmeter 
was  used  to  measure  both  the  voltage  across  the  resistor  and 
the  voltage  across  the  cartridge  heater.  All  the  thermo¬ 
couples  were  connected  to  a  computer  data  acquisition  system, 
and  ambient  conditions  and  voltage  readings  were  supplied 
manually.  The  microcomputer  processed  all  this  information 
through  a  program  included  in  Appendix  A. 

D.  CALIBRATION  OF  THERMOCOUPLES 

A  total  of  fifteen  wall  thermocouples  in  tubes  1  and  7 
(one  thermocouple  in  tube  1  was  defective,  hence  it  was 
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AUXILIARY  HEATER 


Figure  1 . 


Schematic  of  Test  Apparatus 


Figure  1  shows  a  schematic  of  the  test  apparatus.  The 
boiling  surfaces  consist  of  a  cylindrical  copper  cube  with 
a  smooth  interior  surface  and  a  commercially-available  T- 
shaped  finned  surface  manufactured  by  Wieland-Werke ,  AG.  Each 
tube  was  conditioned  for  the  insertion  of  thermocouples  by 
drilling  eight  longitudinal  holes  as  indicated  in  Figure  2 
and  introducing  the  respective  copper-constantan  thermocouples 
centered  by  means  of  capillary  copper  tubes.  In  the  center 
of  the  tube  a  cartridge  heater  was  placed  to  act  as  the  heat 
source  as  shown  in  Figure  3. 

Both  ends  of  the  test  section  were  machined  to  remove 
the  enhanced  surface  on  the  outside  and  to  increase  the 
interior  diameter.  This  reduced  the  cross  sectional  area 
and  thereby  minimized  the  heat  transfer  from  the  ends  of  the 
test  section.  A  Teflon  plug  was  inserted  into  the  free  end 
for  sealing  and  heat  transfer  reduction  purposes.  The  other 
end  of  the  test  section  was  connected  to  a  long  copper  'ex¬ 
tension  tube  through  a  90-degree  elbow,  which  allowed  the  test 
section  to  be  placed  in  a  horizontal  position  within  a  thick- 
walled  Pyrex  glass  vessel.  Since  one  of  the  requirements 
of  the  apparatus  was  to  maintain  the  liquid  at  saturation 
temperature  and  atmospheric  pressure,  a  plate  heater  and  a 
vent  tube  were  provided  to  satisfy  these  conditions.  In 
addition  avoiding  excessive  loss  of  vapor,  a  water-cooled 
copper  coil  was  placed  in  the  vapor  space  of  the  vessel,  to 
act  as  a  primary  condenser.  The  vessel  was  covered  with  a 
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EXPERIMENTAL  APPARATUS 
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A.  REQUIREMENTS  GOVERNING  DESIGN 

The  purpose  of  this  research  was  the  observation  of 
nucleate  pool  boiling  from  different  Gewa-T  surfaces.  In 
order  to  get  a  more  comprehensive  knowledge  about  the  physi¬ 
cal  mechanisms  that  govern  this  phenomenon,  it  is  necessary 
to  measure  the  parameters  involved  which  affect  the  nucleate 
boiling  process.  The  selected  parameters  to  be  obtained 
are  the  following: 

1.  Barometric  pressure 

2.  Test  section  wall  temperature 

3.  Fluid  bulk  temperature 

4 .  Vapor  temperature 

5.  Surface  heat  flux 

B.  TEST  APPARATUS 

The  experimental  apparatus  designed  by  Lepere  [Ref.  9], 
and  modified  by  Hernandez  [Ref.  8],  with  further  minor  modi¬ 
fications,  was  used  during  this  thesis  work.  To  gather  data 
more  conveniently  at  a  faster  rate,  a  Hewlett  Packard  model 
9826  computer  and  an  HP-3054A  automatic  data-acquis ition 
system  were  added  to  the  experimental  facility.  Lepere 
[Ref.  9]  and  Hernandez  [Ref.  8]  provide  detailed  descriptions 
of  the  apparatus;  therefore,  only  a  brief  description  is 
presented  in  this  thesis. 
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to  improve  this  performance  through  the  use  of  further 
enhancing  techniques. 

There  were  four  major  objectives  in  this  thesis  study: 

1.  Take  data  to  verify  results  of  Hernandez  [Ref.  8]  on 
a  standard  21.2  mm  O.D.  Gewa-T  tube  with  and  without  shroud 

2.  Use  wire  wraps  in  the  inter-fin  area  to  improve 
nucleate-boiling  heat- transfer  performance  from  Gewa-T 
tubes . 

3.  Take  data  using  the  combination  of  best  shroud  (found 
in  1)  and  the  best  number  of  wire  wraps  (found  in  2)  on 
Gewa-T  tubes. 

4.  Take  data  to  study  the  effect  of  fin-tip  spacing  and 
fin  density  on  the  heat- transfer  performance  of  Gewa-T 


tubes . 


increase  in  the  heat-transfer  coefficient  when  compared  to 

plain  surface.  The  High  Flux  surface  was  most  effective  over 

a  broad  range  of  heat  fluxes.  The  Thermoexcel-E  surface 

showed  similar  gains  in  heat- transfer  coefficient  to  that 

2 

of  the  High  Flux  surface  below  10  kW/m  .  The  Gewa-T  surface 

was  not  as  effective  as  the  other  surfaces  at  low  heat  fluxes 

2 

but  performed  comparably  at  high  heat  fluxes  (near  100  kW/m  ) 
Marto  and  Hernandez  [Ref.  7]  studied  the  boiling  per¬ 
formance  of  a  Gewa-T  finned  surface  in  Freon-113.  They 
examined  the  influence  of  placing  aluminum  shrouds  around  the 
test  tubes,  and  found  that  the  thermal  performance  of  these 
surfaces  were  dependent  upon  the  liquid-vapor  flow  within 
the  channels  between  neighboring  T-shaped  fins.  They  also 
examined  the  influence  of  the  fins  by  progressively  machining 
away  the  fin  height  to  arrive  at  a  smooth  cylindrical  surface 
They  reported  that  the  addition  of  the  "T-caps"  to  straight 
fins  produced  the  most  significant  improvement  in  heat- 
transfer  performance  when  compared  to  a  smooth  tube.  The 
use  of  shrouds  increased  the  heat-transf er  coefficient  at 
all  practical  heat  fluxes,  but  decreased  the  heat- transfer 
coefficient  at  very  high  heat  fluxes  when  compared  to  the 
unshrouded  surface. 

B.  THESIS  OBJECTIVES 

Based  on  the  foregoing  discussion,  the  major  motivation 
in  this  work  was  to  further  study  the  boiling  heat- transfer 
performance  of  Gewa-T  surfaces  in  R-113,  and  to  attempt 


development  of  generalized  correlations  or  successful 
theoretical  models. 

Yilmaz,  et  al.  [Ref.  5]  compared  the  nucleate  pool-boiling 
heat-transf er  performance  of  three  copper  tubes:  Linde  High 
Flux,  Wieland  Gewa-T  and  Hitachi  Thermoexcel-E .  Their  data 
taken  with  p-xylene  showed  that  while  all  tubes  performed 
better  than  smooth  tubes,  the  best  performance  was  with  the 
Linde  High  Flux  surface.  The  Thermoexcel-E  tube  performed 
better  than  the  Gewa-T  tube  at  low  heat  flux  and  both  per¬ 
formed  similarly  at  high  heat  flux.  These  data  were  taken 
under  the  fully-established  boiling  regime;  thus,  no  infor¬ 
mation  was  provided  on  the  hysteresis  behavior  or  temperature 
overshoot  of  these  surfaces. 

Bergles  and  Chyu  [Ref.  6]  compared  the  nucleate  pool¬ 
boiling  heat-transf er  performance  of  four  different  Linde 
High  Flux  tubes  to  a  plain  copper  tube  in  distilled  water  and 
in  R-113.  They  used  three  different  methods  of  aging  the  test 
surface  prior  to  collecting  data.  They  reported  significant 
temperature  overhsoots  of  both  plain  and  High  Flux  surfaces 
in  R-113.  They  also  reported  that  these  overshoots  were 
sensitive  to  aging,  initial  subcooling  and  rate  of  power 
increase . 

Marto  and  Lepere  [Ref.  1]  obtained  results  on  the  heat- 
transfer  performance  of  three  heat- transfer  surfaces:  Linde 
High  Flux,  Thermoexcel-E,  Gewa-T,  and  a  plain  copper  sur¬ 
face  in  R-113  and  FC-72.  They  reported  a  two  to  tenfold 
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Nishikawa  and  Ito  [Ref.  4]  recently  discussed  two  methods 
to  augment  nucleate  boiling:  (1)  use  teflon-coated  pits  to 
reduce  wettability  of  boiling  liqud,  and  (2)  manufacture 
surfaces  with  numerous  re-entrant  cavities  that  have  the 
ability  to  trap  vapor  and  keep  the  nucleation  sites  active. 
Due  to  the  high  wettability  of  low  surface-tension  fluids 
(such  as  f luorcarbons) ,  the  use  of  teflon-coated  pits  has  not 
been  successful.  Thus,  the  surfaces  with  re-entrant- type 
cavities  are  the  viable  candidates  for  low  surface-tension 
liquids . 

Taking  advantage  of  the  concept  of  re-entrant  cavities, 
a  number  of  enhanced  boiling  surfaces  have  been  invented 
during  the  last  two  decades.  They  include  the  Linde  High 
Flux  surface  of  Union  Carbide,  Gewa-T  surface  of  Wieland 
and  Thermoexcel-E  surface  of  Hitachi.  All  of  these  surfaces 
are  provided  with  re-entrant  cavities  through  specialized 
manufacturing  processes.  These  re-entrant  cavities  introduce 
yet  more  parameters  to  be  desired. 

In  view  of  the  inadequate  understanding  of  the  nucleate 
pool  boiling  process  on  smooth  tubes,  a  thorough  understand¬ 
ing  of  the  boiling  process  on  advanced  surfaces  may  take 
many  more  years  to  come.  Thus,  pfesent-day  designers  have 
very  little  information  when  considering  advanced  boiling 
surfaces.  It  is  very  important  that  researchers  gather  a 
vast  amount  of  experimental  data  covering  a  wide  range  of 
heat-transf er  fluids.  Such  extensive  data  will  enable  the 


good  indication  of  the  complex  nature  of  the  nucleate  boil¬ 
ing  process.  Chongrungreon  and  Sauer  [Ref.  2]  summarized 
nine(9)  correlations  for  nucleate  pool  boiling  of  single¬ 
component  fluids  from  smooth  tubes.  They  stated  the  inade¬ 
quacy  of  any  of  these  correlations  when  compared  with 
experimental  data. 

Nucleate  pool  boiling  is  strongly  dependent  on  the  sur¬ 
face  structure;  i.e.,  the  heat  flux  increases  with  increase 
in  stable  nucleation  site  density.  The  size,  shape  and  den¬ 
sity  of  such  nucleation  sites  are  almost  impossible  to  pre¬ 
dict.  However,  the  entire  nucleate  pool-boiling  process  is 
dependent  on  the  performance  at  these  nucleation  sites.  It 
is  well  known  that  the  ability  of  these  sites  to  develop  and 
trap  small  vapor  bubbles  is  the  reason  for  improved  performance 
in  nucleate  pool  boiling  in  comparison  with  natural-convective 
boiling. 

When  a  bubble  is  being  generated  at  a  surface  nucleation 
site  (or  a  pore) ,  the  wall  superheat  is  given  by  the  follow¬ 
ing  approximate  theoretical  equation  [Ref.  3]: 

2  o  T 

T  -  T  =  - — 

«  S  ov  hFG  R 

Thus,  for  a  given  fluid,  the  wall  superheat  can  be  de¬ 
creased  (boiler  performance  increases  with  decrease  in  wall 
superheat)  with  increase  in  pore  radius.  Of  course,  this 
radius  must  have  a  practical  upper  limit  since  too  large  a 
pore  would  not  trap  a  vapor  bubble. 


I.  INTRODUCTION 

A.  BACKGROUND 

The  energy  crisis  that  !  -aan  in  1974  has  made  a  consider¬ 
able  impact  on  the  design  and  operation  of  energy-conversion 
devices.  Today,  designers  are  challenged  to  design  more 
and  more  efficient  heat- transfer  systems.  Among  the  numer¬ 
ous  areas  that  designers  are  faced  with  this  challenge, 
boilers  or  evaporators  used  in  many  engineering  applications 
demand  considerable  attention.  The  use  of  advanced  boiling 
surfaces  may  lead  to  significant  improvements  in  efficiency 
and/or  significant  savings  in  the  initial  cost  and  equipment 
sizes.  Also,  due  to  the  ever-decreasing  sizes  of  electronic 
components,  the  presence  of  very  high  heat  fluxes  has  made 
the  cooling  of  such  components  a  challenging  problem.  Marto 
and  Lepere  [Ref.  1]  stated  that  the  use  of  advanced  boiling 
surfaces  may  have  potential  promise  in  the  field  of  electronic 
cooling.  However,  the  lack  of  generalized  performance  data 
and  useful  theoretical  models  for  such  advanced  surfaces 
have  still  kept  designers  away  from  the  utmost  use  of  these 
surfaces . 

Due  to  a  lack  of  complete  understanding,  the  theoretical 
treatment  of  the  nucleate  pool-boiling  process,  even  on  a 
smooth  surface  ,  is  extremely  complex.  The  existence  of  a 
large  number  of  semi-empirical  formulas  is  probably  a  very 
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calibration  line.  Based  on  this  fact,  no  calibration  was 
performed  for  the  other  five  tubes.  Instead,  the  calibration 
curve  generated  for  tubes  1  and  7  was  assumed  for  all  tubes. 


Ill . 


EXPERIMENTAL  PROCEDURES 


A.  NORMAL  OPERATION 

The  Procedure  B  employed  by  Marto  and  Hernandez  [Ref.  7] 
has  been  demonstrated  to  be  the  most  adequate  to  carry  out 
the  boiling  runs.  Because  of  the  large  influence  of  surface 
past  history  on  boiling  incipience,  this  procedure  was  estab¬ 
lished  to  allow  repeatability  and  permit  comparison  of  the 
obtained  data  for  different  enhancements. 

The  procedure  is  as  follows: 

1.  Once  the  test  section  has  been  immersed  in  the  pool, 

it  was  subjected  to  one  hour  of  pre-boiling  by  setting  the 

2 

cartridge  heater  to  give  a  heat  flux  of  about  30  kW/m  .  The 
plate  heater  was  adjusted  to  the  maximum  temperature  for  this 
one-hour  period.  This  initial  vigorous  boiling  served  to 
degas  the  liquid  and  to  force  the  noncondensable  gases  out  to 
the  atmosphere. 

2.  After  the  initial  aging  process,  both  test-section  heater 
and  plate  heater  were  secured,  allowing  cooling  of  the  liquid 
and  the  test  section  for  about  30  minutes. 

3.  The  plate-heater  voltage  was  then  adjusted  in  order  to 
maintain  the  liquid  at  saturated  conditions  and  the  cartridge 
heater  was  turned  on  to  the  first  setting  of  12  volts. 

4.  For  all  the  consecutive  settings  during  both  increasing 
and  decreasing  heat  fluxes,  the  boiling  was  allowed  to 
stabilize  for  five  minutes  at  each  power  setting. 


5.  At  a  given  power  setting,  the  following  data  were 
recorded:  heater  voltage,  precision-resistor  voltage,  vapor 

temperature,  liquid  bulk  temperature  and  wall  temperatures 
of  the  test  section. 

B.  GEWA-T  SURFACES 

•  Three  groups  of  surfaces  were  tested.  The  first  surface 
was  the  same  as  used  by  Hernandez  [Ref.  8]  with  a  fin  density 
of  0.75  fins/mm,  an  outside  diameter  of  21.2  mm,  and  an 
inter-fin  gap  of  0.25  mm.  The  second  group  of  three  tubes 
had  0.75  fins/mm,  an  outside  diameter  of  23.1  mm  and  inter- 
fin  gaps  of  0.15,  0.25  and  0.35,  and  the  third  group  of  three 
tubes  had  1.02  fins/mm  and  all  other  specif ications  were  the 
same  as  the  second  group. 

The  geometry  of  the  enhancement  of  these  seven  Gewa-T 
surfaces  was  observed  by  making  a  longitudinal  cut  of  the 
tubes  and  photographing  them  with  50X  magnification  (Figure  5) 
As  can  be  seen,  the  tubes  with  0.75  fins/mm  had  an  approxi¬ 
mately  circular  cavity  shape,  while  the  1.02  fins/mm  tubes 
had  a  somewhat  elongated  cavity  shape.  This  difference  in 
shape  was  probably  due  to  the  technique  used  in  construction 
to  obtain  the  various  specifications.  All  the  T-shaped  fins 
had  a  fin  height  of  approximately  0.95  mm. 

The  exact  determination  of  diameters  at  the  base  of  the 
fins  was  carried  out  from  the  magnified  photographs  of  the 
cavities.  This  diameter  was  determined  by  subtracting  twice 
the  fin  height  from  the  measured  outside  diameter: 
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TABLE  1 

Dimensions  of  Gewa-T  Finned  Tubes 


TUBE 

NO. 

FINS/ 

mm 

FIN  TIP 

GAP  (mm) 

FIN  HEIGHT 
(mm) 

O.D. 

(mm) 

D2 

(mm) 

1 

0 .75 

0 .25 

0 . 865 

21.2 

19.47 

2 

0.75 

0 .15 

0.980 

24.99 

23.05 

3 

0.75 

0 .25 

0.960 

24.99 

23.07 

4 

0.75 

0 . 35 

0.970 

24.99 

23.05 

5 

1.02 

0 . 15 

1.050 

25.30 

23.20 

6 

1.02 

0.25 

1.050 

25.30 

23.20 

7 

1.02 

0.35 

1.125 

25.30 

23.05 

For  simplicity,  the  diameter  at  the  base  (D£)  was  assumed 
to  be  19.5  mm  for  Tube  1  and  23.1  mm  for  all  other  tubes. 

C.  HEAT-FLUX  CALCULATION 

According  to  earlier  description  about  test-section  con¬ 
struction,  the  cartridge  heater  acts  as  the  heat  source, 
inserted  into  the  tube,  and  variations  in  heat- transfer  rate 
are  made  through  different  power  settings  with  a  variac. 

Q„  =  VI  (1) 

H 
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Q  =  heat-transfer  rate  from  the  cartridge 
heater  (W) 

V  =  voltage  across  the  cartridge  heater 
element  (volts) 

I  =  current  through  the  heater  element  (amps) 

This  circuit  is  connected  in  series  with  a  precision 
resistor,  which  provides  the  means  to  compute  the  electrical 
current  by  the  measurement  of  the  voltage  drop  across  the 
res istor : 

1  -  vres/e  (2> 

where : 

Vres  =  voltage  drop  across  the  precision 
resistor  (volts) 

R  =  resistance  of  the  precision  resistor  (ohms) 

=  2.031  ohms 

Both  ends  of  the  Gewa-T  tube  were  machined  to  provide  a 
smooth  surface  and  a  reduced  wall  thickness  in  order  to  mini¬ 
mize  longitudinal  heat  conduction  and  therefore  the  heat 
losses  from  the  unenhanced  ends.  The  geometry  of  these  ends 
is  shown  in  Figure  6,  where: 

b  =  base  of  the  T-shaped  fins 

T,  =  temperature  at  the  base  of  T-shaped  fins 


where  : 

2 

Ac  =  cross-sectional  area  of  tube  wall  (m  ) 

Dq  =  tube  outside  diameter  (m) 

=  tube  inside  diameter  (m) 

2 

Ag  =  tube  outside  surface  area  (m  ) 

P  =  tube  outside  wall  perimeter  (m) 

It  was  assumed  that  the  temperature  at  the  base  of  the 
thin-walled  fin  was  equal  to  the  temperature  at  the  base  of 
the  T-shaped  fins.  Therefore: 
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e  (x) 


T  (x)  -  T 


SAT 


(5) 


In  general, 


d29  2. 

— ~  -mo 

dx 


=  0 


(6) 


where : 


m 


h  P 
kA„ 


(7) 


Therefore 


QF  qb 


=  -kA 


dT 
C  dx 


=  -kA 


x=0 


d8  ! 

C  dx 


x=0 


(s; 


Now,  assuming  negligible  convective  heat  loss  from  the 
thin  tube  tip: 


d0 


0 


(9) 


and 


Cosh  m(L-x) 
Cosh  mL 


(10) 


For  this  case 


Qf  =  4  Pkc  Ac  eB  Tanh  ml^ 


(ID 


34 


(o) 


=  T  -  T 

B  SAT 


(12) 


where 


'SAT 


=  heat-transfer  rate  through  unenhanced 
end  (W) 

=  fluid  saturation  temperature  (°C) 

=  thermal  conductivity  of  the  copper  (W/m.K) 

2 

=  average  heat- transfer  coefficient  (W/m  .K) 


Assuming  the  wall  temperature  is  the  same  as  the  average 
of  the  eight  calibrated  wall  temperature  thermocouple 
measurements,  less  temperature  drop  due  to  radial  conduction 
from  thermocouple  position  to  the  base  of  the  T-shaped  fin, 


1 n(D2/Dx) 


TAVG  QH  2  7T  L  k” 


13) 


AVG 


8 

(  l  T  )/8.0 
n=l  n 


(14) 


where 


D1  =  diameter  at  the  position  of  thermocouples 
1  (m) 

T  =  average  wall  temperature  (°C) 

AV  Cj 


For  the  geometry  shown  in  Figure  6: 


T  t 

LU  +  2 
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(15) 


t  =  D  -  D. 

o  1 


(16) 


where  : 

=  corrected  length  (m) 

=  length  of  unenhanced  ends  (m) 
t  =  tube  thickness  (m) 

The  average  difference  between  the  wall  temperature  and 
the  saturation  temperature  may  now  be  determined  from  the 
following  equation: 


,  C  Cosh  m(L--x) 

=  7 L  J  C 


C  0 


Cosh  mL, 


dx 


(17) 


0 


mLC 


Tanh 


mLc 


(18) 


0  = 


T  -  T 
W  SAT 


(  19) 


where : 


(TW~T  T)  =  average  difference  between  surface 
'  wall  temperature  and  fluid 

saturation  temperature  (K) 


Since  the  unenhanced  ends  are  very  long  and  the  thickness 
is  very  small,  only  very  little  heat  transfer  would  take 
place  at  these  ends.  Thus,  these  ends  undergo  free  convection, 
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for  which  the  Churchill  and  Chu  correlation  as  stated  in 


'Ref.  10]  for  the  average  Nusselt  number  was  assumed: 


Nur 


0.387  Ra 


1/6 

D 


0.60+- 


;i  +  (0.559/Pr) 


9/16  j 8/27 


-5  -  12 

10  <  RaD  <  10 

o 


The  average  Nusselt  number  is 


h  D 


Nu 


D 


where : 


k  =  thermal  conductivity  of  Freon-113  (W/m.K) 


Solving  for  h: 


*  0.387  RaJ/6 

H  =  t  |°'60+7l  +  (0.559/Pr)9/16]8/27 


Pr 


v_ 

a 


Rar 


q  3  (T  -  T  )  D 
y  v  W  SAT  o 


v  i 


1  up 

IT 
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22) 


(23) 


(24) 


(25) 


.  .  JLj»  *  m.  f. 


where  : 


Pr 


Ra. 


Prandtl  number 

2 

kinematic  viscosity  (m  /s) 

2 

thermal  diffusivity  (m  /s) 
average  Rayleigh  number 

2 

gravitational  acceleration  (m/s  ) 

volumetric  thermal  expansion  coefficient 
( 1/K) 

density  of  freon-113  (kg/m^) 


Now,  substitution  of  Equations  7,  18  and  24  into  Equation  22 
results  in: 


0.387 


1/2 


1.1 


JCl 


'  2LC(kA^}  1/2 


1/6 


(26) 


h  =  ;o.60  + 

o 


[1  +  (0.559/Pr) 


9/16, 8/27 


Equation  26  is  solved  for  h  by  iterative  technique 
within  a  range  of  precision  of  0.001.  Knowing  the  value  of 
natural-convective  heat- trans f er  coefficient  at  unenhanced 
ends,  the  total  heat-loss  rate  is  calculated  from  Equation 
11,  and  the  hea t- t r ans f er  rate  through  the  enhanced  Gewa-T 
surface  is  obtained  by  subtracting  the  total  heat- loss  rate 
from  the  total  heat-transfer  rate: 


O. 


'  LOSS 


=  20. 


(27) 


38 


diameter.  For  this  purpose,  a  semi-empir xcal  correlation 
develooed  bv  Cornwell  et  al.  [Ref.  11]  is  used: 


CTB  Re™3 


hTB  D 


ID  .  2/3 


'TB  hFG  V 


D  V  A 


G  =  q/h 


CTB  kF,hTB  1TD,2-/3 
0  kFG “F 


CTB  kF  ,  AT  D  2/3 
D  hFG  ^F 


(^)3  (^^)2 
hFG  J  F 


(CTB  kF)3  (hFG\;F!  D 


Using  this  relationship,  Hernandez  *s  data  were  revised 
and  a  least-squares-f it  curve  was  generated.  This  curve  is 
also  plotted  in  Figure  12. 


H.  BUILING  PERFORMANCE  OF  GEWA-T  TUBE  WITH  SHROUDS 


The  same  four  shrouds  tested  by  Hernandez  were  used 
during  this  investigation.  These  shrouds  had  the  following 


geometry : 


agree  to  within  0.2  percent  at  a  heat  flux  of  30,000  to 
2 

100,000  W/m  .  These  data  runs  show  considerably  good  repeata¬ 
bility  of  this  experiment.  (Similar  agreements  were  also 
found  for  all  other  tubes.) 

Also,  these  data  show  a  considerable  enhancement  of 
Gewa-T  tube  performance  over  the  "Smooth"  tube  data  of 
Hernandez.  For  example,  at  IT  =  3.1  K,  the  Gewa-T  tube 
gives  a  4.65  times  larger  heat  flux  than  the  smooth  tube. 

As  also  discussed  by  Lepere  [Ref.  9]  and  Hernandez  [Ref. 

8]  this  enhancement  in  boiling  performance  is  mainly  attri¬ 
buted  to  the  existence  of  re-entrant  cavities  in  the  Gewa-T 
tube . 

G.  BOILING  PERFORMANCE  OF  SMOOTH  TUBE 

To  compare  the  performance  of  Gewa-T  tubes  to  smooth-tube 
performance,  data  of  Hernandez  were  chosen.  It  is  important 
to  point  out  that  his  smooth  tube  did  not  consist  of  a  commer¬ 
cially-  available  surface.  Instead,  Hernandez  obtained  his 
"Smooth"  surface  by  machining  fins  off  of  a  Gewa-T  tube. 

He  did  not  remove  fins  completely;  the  remaining  portions  on 
the  tube  measured  about  0.1  mm. 

Figure  12  shows  Hernandez 's data  (heat  flux  as  a  function 
of  wall  superheat)  for  only  decreasing  heat  fluxes.  A 
least-squares- f it  curve  generated  for  these  data  points  will 
be  used  in  future  figures  representing  data  on  tube  1.  Since 
tubes  2  through  7  have  a  larger  diameter,  it  was  necessary 
to  revise  Hernandez's  data  to  include  the  difference  in 
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Figure  11.  Reproducibility  Test:  Plain  Gewa-T  Tube 


E.  REPEATABILITY  WITH  PREVIOUS  DATA 

Following  modifications  made,  especially  to  the  data 
gathering  means  (i.e.,  automatic  data-acquisition  system  and 
newly-written  computer  programs) ,  it  was  necessary  to  take 
data  for  tube  1  (this  is  the  same  tube  tested  by  Hernandez) 
to  check  the  present  data  with  Hernandez  data  under  similar 
conditions.  Even  though  a  total  of  5  different  runs  were 
made  during  this  stage  through  debugging  of  programs  and 
system  familiarization,  only  run  number  5  will  be  presented 
in  this  thesis  with  a  comparison  to  Hernandez  data. 

Figure  10  shows  the  comparison  between  the  present  data 
and  the  data  of  Hernandez.  A  very  good  agreement  is  evident 
as  the  two  sets  disagree  by  less  than  2  percent  for  the 
nucleate  boiling  regime.  However,  considerable  disagreement 
exists  in  the  region  of  increasing  heat  flux  before  total 
nucleate  boiling.  This  is  believed  to  be  due  to  the  rather 
random  behavior  of  the  boiling  process  during  this  regime 
as  it  will  be  presented  in  the  next  section,  even  two  consecu¬ 
tive  runs  for  a  given  tube  would  not  give  the  same  results 
during  the  unstable  boiling  regime. 

F.  REPEATABILITY  STUDY  OF  PRESENT  WORK 

Figure  11  shows  data  from  two  different  runs  completed 
on  two  different  days.  As  discussed  earlier,  a  considerable 
disagreement  can  be  seen  before  the  tube  undergoes  complete 
nucleate-boiling  process.  However,  for  the  complete 
nucleate-boiling  process,  the  data  from  different  data  runs 
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confirmed  the  validity  of  the  mathematical  model  assumed  for 
data  reduction.  However,  a  small  rate  of  bubble  formation 
was  evident  at  the  tip  of  the  unenhanced  fin  end;  these 
bubbles  were  in  fact  generated  at  the  base  of  the  fin  (i.e., 
between  the  Teflon  plug  and  the  inner  surface  of  the  fin) . 

The  vapor  bubbles  generated  in  the  inter-fin  cavities 
are  carried  up  into  the  channels  around  the  tube  due  to 
buoyancy  forces.  As  these  bubbles  travel  up,  they  activate 
more  and  more  nucleation  sites  as  well  as  they  help  for  the 
removal  of  larger  bubbles  from  the  channels.  This  phenomenon 
is  generally  referred  to  as  the  “Chimney  Effect." 

If  the  heat  flux  is  too  large  (see  Figure  7),  the  rate  of 
vapor  bubble  generation  may  be  too  large  that  the  interfin 
channels  could  be  full  of  a  vapor  blanket  that  may  disrupt 
the  flow  of  liquid  to  the  hot  metal  surface  where  active 
boiling  takes  places.  At  this  point,  further  increase  in 
heat  flux  may  result  in  excessive  wall  temperatures.  Under 
these  conditions,  vapor  bubbles  were  seen  ejecting  radially 
outward  from  the  test  surface.  Line  D-E  represents  the  boil¬ 
ing  regime  just  described.  It  is  clear  that  the  slope  (i.e., 
heat- transfer  coefficient)  at  E  is  considerably  smaller  than 
that  at  point  D.  This  decrease  in  heat-transf er  coefficient 
is  assumed  to  be  due  to  the  existence  of  an  insulating  vapor 
blanket . 

When  heat  flux  is  gradually  decreased,  the  curve  follows 
a  different  path  (D-F)  as  shown  in  Figure  7.  This  is  due 
to  the  existence  of  stable  nucleation  sites  that  will  remain 
active  for  a  wide  range  of  heat  fluxes. 
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yure  9.  Photograph  of  Gewa-T  Surface 


Figure  8.  Photograph  of  Gewa-T  Surface  at  High  Ilea 


From  point  A  to  point  B,  a  continuous  increase  in  heat 


flux,  accompanied  by  an  increase  in  wall  superheat  (Tw  -  TSAT) 
is  observed  when  heat  flux  is  increased.  This  region  corres¬ 
ponds  to  natural-convection  process,  as  also  verified  during 
experimental  runs  with  the  absence  of  bubbles  along  both 
enhanced  surface  and  unenhanced  ends. 

From  point  B  to  point  C,  a  dfastic  reduction  in  wall 
superheat  was  observed  while  the  heat  flux  continuously 
increased.  This  region  was  somewhat  unstable  as  more  and 
more  nucleation  sites  became  active.  This  region  is  known 
as  the  mixed  boiling  region,  where  transition  from  natural- 
convective  boiling  to  nucleate  pool  boiling  takes  place.  At 
these  heat  fluxes,  visual  observations  revealed  that  the 
bubble  formation  was  quite  random.  While  some  portions  of 
the  tube  showed  bubbles,  other  portions  showed  no  bubble 
formation.  The  wall  superheat  continued  to  decrease  until 
all  nucleation  sites  became  active,  as  indicated  by  point  D. 

Once  all  available  nucleation  sites  have  been  activated, 
the  wall  superheat  starts  to  increase  with  increasing  heat 
flux  as  shown  by  region  D  to  E.  When  boiling  in  this  region, 
the  entire  tube  showed  very  high  density  of  bubble  formation. 

At  this  point,  it  is  important  to  point  out  that  during  all 
experimental  runs,  for  the  entire  range  of  power  settings  (7.5  to 
4  8  5  W)  ,  little  or  no  nucleate  pool  boiling  on  the  unenhanced 
surface  of  the  ends  was  observed  (see  Figures  8  and  9)  for 
the  different  test  sections  tested.  This  observation 
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AT 


A-B:  Natural  Convection 

B-C:  Mixed  Boiling  (Activation  of 

Nucleation  Sites) 

C-D:  Active  Boiling 

D-E:  High  Flux  Nucleate  Boiling 

E-F:  Stable  Nucleate  Boiling 


Figure  7. 


Typical  Nucleate  Pool-Boiling  Curve 
for  a  Gewa-T  Surface 


runs  were  made.  All  of  these  runs  showed  only  a  negligible 
variation  of  the  wall  thermocouple  temperatures;  maximum 
variation  was  t0.25  K.  These  runs  clearly  showed  that  the 
heater  element  indeed  provided  a  uniform  heat  flux  along  the 
circumferential  direction.  No  detailed  data  are  presented 
in  this  thesis  for  the  tube  in  the  vertical  position.  All 
other  results  presented  herein  are  for  the  tubes  in  the 
horizontal  position. 

C.  LONGITUDINAL  UNIFORMITY  OF  HEAT  FLUX  ALONG  THE  TEST 

SECTION 

During  most  of  the  data  runs  with  the  tubes  in  the  hori¬ 
zontal  position,  somewhat  unexplainable  phenomena  were 
observed.  As  the  heat  flux  was  gradually  increased,  bubble 
formation  began  simultaneously  at  the  two  ends,  while  the 
center  portion  of  the  tube  underwent  only  natural  convection. 
Approximately  40  percent  of  maximum  heat  flux  was  necessary 
to  obtain  nucleate  boiling  from  the  entire  tube.  Since  all 
the  wall  thermocouples  are  located  at  midway  along  the  tube, 
the  average  wall  temperature  values  provided  in  all  data 
runs  may  be  in  slight  error. 

D.  PLOT  ANALYSIS  FOR  PLAIN  SURFACE 

Figure  7  represents  a  typical  curve  for  nucleate  pool¬ 
boiling  process  of  Gewa-T  surface  in  Freon-113.  The  observed 
behavior  of  this  process  is  analyzed  in  this  section  by 
studying  the  different  sections  of  the  boiling  curve. 
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separate  file.  This  plot  data  file  could  easily  be  recognized 
by  the  "P"  added  at  the  beginning  of  the  raw  data  file  name. 

In  this  manner,  the  plot  data  file  corresponding  to  the  above 
raw  data  file  has  the  name  PGT5H28. 

A  total  of  thirty  different  data  runs  were  carried  out 
during  this  investigation.  These  runs  break  down  in  the 
following  manner: 


Run  Number 

Tube  Number 

1-3 

1 

4-10 

1 

11  -  21 

1 

22  -  23 

2 

24  -  25 

3 

26 

4 

27  -  28 

5 

29 

6 

30 

7 

Test  Condition 
Plain  -  Vertical  position 
Plain  -  Horizontal 
With  shrouds  and/or  Wires 
Plain  -  Horizontal 
Plain  -  Horizontal 
Plain  -  Horizontal 
Plain  -  Horizontal 
Plain  -  Horizontal 
Plain  -  Horizontal 


B.  CIRCUMFERENTIAL  UNIFORMITY  OF  HEAT  FLUX  AROUND  HEATER 
ELEMENT 

During  the  early  stages  of  this  investigation,  it  was 
thought  necessary  to  test  if  the  heater  element  would  provide 
a  non-uniform  heat  flux  along  the  circumferential  direction. 
Such  a  non-uniform  heat-flux  variation  was  assumed  possible 
based  on  the  heater-element  construction.  Therefore,  tube  1 
was  immersed  vertically  in  the  R-113  liquid  pool  and  three 
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IV.  RESULTS  AND  DISCUSSION 

A.  DESIGNATION  OF  DATA  FILES 

A  total  of  seven  Gewa-T  tubes  were  tested  during  this 
investigation.  As  can  be  seen  from  Table  1,  the  primary 
geometric  variables  were:  fin-tip  gap,  fin  density  and 
outside  diameter  of  tubes  including  fins. 

The  general  designation  for  data  files  is: 

GTNPnn 


where : 


GT  indicates  Gewa-T  surface 

N  specifies  the  number  of  the  tube  tested 
(Table  1) 

P  indicates  the  position  of  the  test  section 

(i.e.,  H  for  horizontal  and  V  for  vertical), 
and 

nn  indicates  the  corresponding  cumulative  run 
number . 

For  example,  GT5H28  specifies  tube  number  five  tested  in  the 
horizontal  position  during  run  number  28.  From  Table  1, 
the  corresponding  specifications  for  this  test  section  are: 
25.3  mm  outside  diameter,  23.1  mm  diameter  at  the  base  of  the 
fins,  1.02  fins/mm  and  inter-fin  gap  of  0.15  mm.  While  the 
file  name  just  presented  contain  data  in  the  very  raw  form, 
the  processed  data  (i.e.,  q  versus  AT)  were  stored  in  a 


40 


Shroud  Number 
1 
2 

3 

4 


Upper  Opening 
60° 

60° 

30° 

60° 


Lower  Opening 
60° 

30  0 
30  3 
8  .  5° 


All  of  these  shrouds  had  been  designed  for  tube  1;  thus, 
no  shrouds  were  used  for  tubes  2  through  7 . 

Figures  13,  14,  15,  and  16  show  data  on  tube  1  with  shroud 
numbers  1,  2,  3  and  4,  respectively.  As  also  reported  by 
Hernandez  [Ref.  8],  all  the  shrouds  improved  boiling  per¬ 
formance  at  low  heat  fluxes  compared  to  plain  Gewa-T  tube. 
Table  2  lists  performance  calculations  for  tube  1  with  various 
shrouds.  It  can  be  seen  that  the  60°  ’<■8.5°  shroud  (number  4) 

produced  the  maximum  improvement  of  253  percent  in  the  boiling 

2 

coefficient  over  smooth  tube  at  a  heat  flux  of  10,000  W/m  . 

At  this  heat  flux,  the  30°  *30°  shroud  resulted  in  a  165 

percent  improvement  in  the  bciling  heat-transf er  coefficient. 

2 

At  a  heat  flux  of  50,000  W/m  ,  all  shrouds  resulted  in  about 
190  percent  improvement,  except  the  30°  *30°  shroud  which 
resulted  in  a  180  percent  improvement. 

2 

At  very  high  heat  fluxes,  for  example  at  q  =  100,000  W/m  , 
the  60°  x  60°  shroud  showed  a  maximum  performance  increase  of 
106  percent;  the  60°  x  30°  shroud  resulted  in  a  49  percent 
improvement;  the  60°  <8.5°  shroud  resulted  in  a  18  percent 
increase;  and  the  30°  *30°  shroud  resulted  in  an  13  percent 
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Effect  of  60°  >60°  Shroud  on  Boiling  Performance 
of  Plain  Gewa-T  Tube  1 


ffect  of  60°  x  30°  Shroud  on  Boiling  Performance 


Summary  of  Results 
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improvement.  It  can  be  seen  that  the  values  just  presented 
are  somewhat  proportional  to  the  summation  of  the  two  shroud 
angles.  This  behavior  indicates  that  at  very  high  heat 
fluxes,  too  small  a  shroud  opening  can  obstruct  the  outflow 
of  vapor  and/or  the  inflow  of  liquid. 

However,  it  must  be  noted  that  at  all  practical  heat 

2 

fluxes  (10,000  to  80,000  W/m  ),  all  the  shrouds  resulted  in 
a  considerable  enhancement  in  the  boiling  performance.  The 
60°  x  8.5°  shroud  appears  to  be  the  best  selection  for  all 
practical  heat  fluxes. 

I.  PERFORMANCE  OF  GEWA-T  TUBE  WITH  WIRE  WRAPS  IN  INTER-FIN  GAP 

A  series  of  runs  were  made  with  two  to  five  0.1  mm-diameter, 
copper  wire  wraps  in  the  inter-fin  gap  of  tube  1.  Figures 
17,  18,  19  and  20  show  data  with  wire  wraps  of  2,  3,  4  and 
5,  respectively.  It  can  be  seen  that  all  of  these  curves 
show  increases  in  boiling  performance  at  all  wall  superheats 
for  decreasing  heat  flux  values.  This  enhancement  is  most 
probably  due  to  the  increase  in  the  nucleation  sites  within 
the  inter-fin  cavities.  However,  there  is  a  small  reduction 
in  performance  in  the  natural  convection  process  probably 
owing  to  that  wires  could  act  as  an  insulating  material  due 
to  the  thermal  contact  resistance  between  the  wires  and  the 
inter-fin  area. 

Table  2  also  shows  that  the  use  of  three  wires  resulted 
in  the  best  performance  increase.  This  combination  resulted 
in  a  341  percent  improvement  in  the  boiling  heat- transfer 


Wrapped  in  the  In 
liling  Performance 


Figure  19.  Effect  of  Four  Wires  Wrapped  in  the  Inter-fin 
Cavity  of  Tube  1  on  Boiling  Performance 


coefficient  over  the  smooth  tube  value  at  a  wall  superheat 
of  1.4  K. 

J.  BOILING  PERFORMANCE  OF  TUBE  1  WITH  60°  *8.5°  SHROUD  AND 

THREE  WIRES 

Data  were  also  taken  to  combine  the  best  shroud  and  the 
best  number  of  wire  wraps  in  the  inter-fin  channel.  For 
this  purpose,  three  wire  wraps  were  added  to  the  inter-fin 
channel  and  data  were  taken  with  the  60°  x  8 . 5°  shroud.  As 
Figure  21  and  Table  2  show,  this  combination  resulted  in  a 
performance  yet  superior  to  the  use  of  three  wires  or  the  use 
of  60°  x  8.5°  shroud  separately.  Table  2  shows  that  this 
arrangement  resulted  in  481  percent  increase  in  the  boiling 
heat-transf er  coefficient  over  the  smooth-tube  value  at  a 
heat  flux  of  10,000  W/m2 . 

K.  EFFECTS  OF  FIN-TIP  GAP  AND  FIN  DENSITY  OF  GEWA-T 

TUBES  ON  BOILING  PERFORMANCE 

As  described  earlier,  six  Gewa-T  surfaces  were  tested 
during  this  stage  of  the  investigation;  tubes  2,  3,  and  4 
had  a  fin  density  of  0.75  fins/mm,  and  tubes  5,  6  and  7 
with  a  fin  density  of  1.02  fins/mm,  and  each  set  of  tubes 
with  fin-tip  gaps  of  0.15,  0.25  and  0.33  (see  Table  1). 

Micrographs  were  obtained  for  all  enhanced  surfaces  on 
a  Scanning  Electron  Microscope  with  a  magnification  of  50X. 
Figure  5  shows  micrographs  representing  tubes  2  and  5. 
Computed  surface  areas  and  area  ratios  (with  respect  to 
smooth  tube)  are  listed  in  Table  3. 
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TABLE  3 


Computed  Surface  Areas  and  Area  Ratios 


Tube 


1 

2 

3 

4 

5 

6 
7 


Cavity  and  Fin2 
Surface  Area  (m  ) 

240 .47  x  lo-6 

288.01  x 10~6 

305.34  x lo'6 
285.07  x  io~6 
256.04  x lo“6 
248.61  < 10~6 
267.91  x  10~6 


Total  Surface  Area  Ap/ ASM00^H 

(m2) 


8 . 818  x  io  J  - 

10.944  x  10-3  2.74 

11. 603  x io-3  2.91 

10.833  •  IO"3  2 . 72 

13.314  •  10"3  3.60 

12.928  <  10-3  3.50 

13.931  x  IQ-3  3.77 


It  can  be  seen  that  tubes  with  a  higher  fin  density 
(1.02  fins/mm)  result  in  larger  area  ratios  than  the  tubes 
with  lower  fin  density  (0.75  fins/mm).  However,  for  a  given’ 
fin  density,  the  variation  of  the  area  ratio  with  fin-tip 
gap  is  not  conclusive. 

Figures  22  through  27  show  data  taken  on  tubes  2  through 
7,  respectively.  These  results  are  also  listed  in  Table  2. 
This  table  shows  that  tube  6  (with  a  fin  density  of  1.02 
fins/mm  and  a  fin-tip  gap  of  0.25  mm)  resulted  in  the  highest 
heat  flux  at  all  values  of  wall  superheat. 

To  investigate  the  effects  of  fin  density  and  fin-tip 
gap  on  boiling  heat-transfer  performance,  the  results  pre¬ 
sented  in  Figures  22  through  27  (which  showed  variations  of 


67 


T  i.  i  ■  ■ 


APPENDIX  A 

DATA  REDUCTION  PROGRAM 


,  • '  1 

:  IF 

D"' 

•  |  *  ;‘i  » 

'JftrE  :  Mdv  ’  5  .  * 

■  ■  J  Q  t 

'EVi'iEJ:  j\± 

*  /  JO  .  *  9Ha 

'  ;  41) 

BEE3 

*  ')50 

3R  1  f\j  7  R  T  5  * 

jhO 

3RIN'  .8 IMG  " 

4 X .  Select  option: 

i  1170 

PRINT  USING  " 

SX  .  “'*0  Taxing  data  or  re~o  recess  i  no  previous 

■  ;30 

PRINT  USING  " 

6X  .  '  Plotting  data 

•  090 

I N P 1  j T  loo 

:  '  oo 

IF  I dp  -  0  rHEN 

CALL  Main 

’  1  10 

IF  Ido-1  THEN 

CAL_  Plot 

!  120 

END 

1  130 

SUB  Mam 

1  1  40 

COM  /Co/  <2  C  7  > 

’150 

DIM  Emf < 9 ) .  T ( 9 ) 

'  ISO 

3978  i).  10086091  .  25727  .  .14369  76  7345 . 3295 . 79025595 . 9 1 

'  ’  70 

DATS  -924  7486589.  b.9768i3E*-i  ’  . -2 . 661 92E"  3 . 3 . 940  73E  + 1  4 

!  '90 

READ  C ( *  > 

"90 

BEEP 

1200 

INPUT  •'SELECT 

T1JBE  QD  (0-0  75  in.  '-1  .0  in)”. ltd 

'210 

IF  ltd-0  rHEN 

1  220 

31 * .0159 

’  Diameter  at  thermocouple  positions 

1230 

D2- .0133 

!  Diameter  of  test  section  to  tne  oase  of  *  i n 

'240 

D i * . 0 1 93 

*  .nsice  diameter  of  •mennancod  enas 

1250 

Do-. 01 30 

'  Outside  diameter  o'  unem-ianceo  .-nos 

'250 

L* .0490 

!  Lenqtn  of  ennanced  surface 

'270 

lu= .0319 

!  Length  of  unenhanced  surface  a  the  ends. 

'290 

ELSE 

’  290 

Di - . 01 796 

'  300 

D2- .0231 

'310 

D i - . 02032 

I 

Do- . 02222 

1  330 

l* .0508 

1  140 

L Li-  .031  75 

'350 

end  :r 

•  160 

■_c*Lu*<  Do-0  i  J 

/2 

’  370 

M-o T »<  Dr.  2-Di 

2 )  /4 

i  380 

4-=>r-flr> 

'  390 

<cu=  195 

'  Thermal  conductivity  of  Ioddot 

'400 

°PINTE-  '0 

'4  10 

CLEAR  '09 

’  420 

BEE3 

•430 

INPUT  "ENTER 

MONTH .  DM  IE  AND  TIME  ( MM : DD : HH : MM : SS ) " . Dat eS 

1  441- 

OUT3UT  ’0rj :  ”  T 

D"  :  Date'r. 

’  450 

jlj;p*jr  toq:,,t 

0” 

*  4hfj 

ENTER  '09  :  Oat*??. 

'  j70 

3R  ['1  r  •• 

yontn,  iace  and  time  :":Dat:e5 

•  480 

3R  IN ' 

■  JHQ 

DRINT  USING  " 

0X.  MJ  : L :  3roqram  name  :  BOIL 

■  son 

3EE3 

'  r:  ;  1] 

INP'.i  "  "r‘l 

DISK  ‘lUMBER"  .  Dn 

'520 

PRINT  js;NG  ■ 

1  bX  .  D  l Sk  niinoer  ~  ""  .Z?"  ;C’n 

’  530 

BEE3 

'5411 

INPUT  "E'i rER 

INPUT  mqde  (  o-  30 54P  .  1  IlE  )"  .  I.k 

*  5  0  'J 

».  n  =  .  n  *  * 

“  f,  i  i 

.  n  -  ’  yl 

•  s ; 

• 0) 

INPL1'  "OIUF  4 

NAME  "R  '-‘E  RAW  MATA  3Iu-".D  -  I-S 

•  59  n 

33  IN '  USING  ' 

*  5X  .  t  .  e  .  *  4fl”  ;  J  •  .  .  ?. 

i 


81 


VI .  RECOMMENDATIONS 


(1)  Investigate  the  reason  for  the  existence  of  a 
longitudinal  temperature  profile  and  take  appropriate 
steps  to  alleviate  this  problem. 

(2)  Study  the  possibility  of  adding  roughness  in  the 
inter-fin  cavity  of  Gewa-T  tubes.  This  operation  should 
be  performed  before  the  rolling  process  to  obtain  the  T- 
shaped  fins. 

(3)  Repeat  runs  made  during  this  investigation  with 
water  as  the  boiling  medium. 


This  combination  resulted  in  a  481  percent  increase  in  boil¬ 
ing  heat-transf er  coefficient  at  a  wall  superheat  of  1.1  K. 

(6)  The  use  of  wires  appears  to  be  superior  to  the  use 
of  more  expensive  shrouds. 

(7)  Among  the  fin-tip  gaps  of  0.15,  0.25  and  0.35  mm 
tested,  the  0.25  mm  gap  resulted  in  the  best  boiling 
performance . 

(8)  The  Gewa-T  tubes  with  a  fin  density  of  1.02  fins/mm 
performed  superior  to  a  fin  density  of  0.75  fins/mm.  At  a 
wall  superheat  of  3.5  K,  the  lower-fin-density  tube  increased 
the  boiling  heat-transf er  coefficient  by  80  percent,  while 
the  higher-fin-density  tube  increased  by  103  percent. 


V.  CONCLUSIONS 


(1)  Experimental  observations  revealed  a  non-uniform 
temperature  distribution  along  test  section.  The  bubble 
generation  was  seen  to  start  at  the  two  ends  and  migrate  into 
the  center  as  the  heat  flux  was  increased.  No  clear  explana¬ 
tion  is  possible  for  this  temperature  distribution. 

(2)  Results  of  the  present  investigation  showed  very- 
good  repeatability  with  the  results  of  Hernandez  obtained 
under  similar  conditions. 

(3)  The  use  of  shrouds  improved  the  boiling  heat- transfer 

coefficient  of  Gewa-T  surface  at  all  practical  heat  fluxes. 

However,  a  considerable  reduction  in  heat-transf er  coefficient 

2 

was  seen  at  very  high  heat  fluxes  (>  100,000  W/m  ) .  Among 
the  four  shrouds  tested,  the  shroud  with  60°  *8.5°  openings 
resulted  in  the  best  performance  verifying  Hernandez's 
results.  This  shroud  resulted  in  a  253  percent  increase  in 
boiling  heat-trans fer  coefficient,  over  smooth  tube,  at  a 
wall  superheat  of  1.8  K. 

(4)  The  use  of  0.1  mm-diameter  wires  (2  to  5)  improved 
boiling  heat- transfer  coefficient  at  all  practical  heat 
fluxes.  The  use  of  3  wires  resulted  in  the  best  performance 
with  a  341  percent  increase  in  the  heat- transfer  coefficient 
(over  the  smooth-tube  value)  at  a  wall  superheat  of  1.4  K. 

(5)  Further  improvement  in  boiling  performance  was 
obtained  by  using  the  60°  x8.5°  shroud  with  three  wire  wraps. 
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explanations  given  above  should  be  valid  for  the  poorer  per¬ 
formance  of  cubes  2  and  4. 

It  can  be  seen  that  all  tubes  show  smaller  enhancement 
ratios  at  low  wall  superheats,  and  the  above-mentioned 
differences  in  performance  are  more  pronounced.  On  the 
other  hand,  all  the  tubes  show  greater  enhancements  at  larger 
wall  superheats  and  they  seem  to  approach  a  common  enhance¬ 
ment  ratio.  The  increase  in  enhancement  is  believed  to  be 
due  to  the  increased  chimney  effect,  which  helps  to  remove 
vapor  bubbles  from  channels  resulting  in  a  thinner  liquid 
film  at  the  hot  boiling  surface.  When  heat  flux  is  further 
increased,  vapor  can  provide  a  blanket  around  the  hot 
surface,  which  may  eventually  result  in  a  critical  heat 
flux.  It  appears  that  the  critical  heat  flux  is  not  strongly 
dependent  on  the  fin  density  or  fin-tip  gap. 

The  results  shown  in  Figure  28  were  also  replotted  as 
shown  in  Figure  29  to  further  study  the  performance  of  these 
six  tubes.  This  figure  shows  the  variation  in  enhancement 
ratio  obtained  beyond  the  area  increase  (i.e.,  the  ratio  of 
heat-transfer  enhancement/area  ratio)  as  a  function  of  wall 
superheat.  It  is  evident  that  obviously  all  tubes  resulted 
in  enhancement  ratios  beyond  the  area  increase. 

In  summary,  a  fin  density  of  1.02  fins/mm  resulted  in 
performance  superior  to  a  fin  density  of  0.75  fins/mm.  Also, 
in  both  cases,  a  fin-tip  gap  of  0.25  mm  resulted  in  the  best 
performance . 
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heat  flux  as  a  function  of  wall  superheat)  were  re-plotted 
on  a  different  basis  as  shown  in  Figure  28.  This  figure  shows 
the  variations  of  the  enhancement  of  heat-transf er  coefficient, 
over  the  smooth-tube  value,  with  wall  superheat  for  all  six 
tubes.  It  is  evident  that  tube  6,  which  has  a  fin  density 
of  1.02  fins/mm  and  a  fin-tip  gap  of  0.25  mm  results  in  the 
maximum  enhancement  at  all  wall  superheat  values. 

Tube  5  shows  poorer  performance  than  tube  6.  This  poorer 
performance  is  believed  to  be  due  to  the  fact  that  the  smaller 
fin-tip  gap  disrupts  the  flow  of  liquid  into  the  channels 
and/or  the  vapor  flow  out  of  the  channels.  Tube  7  (which  has 
the  largest  fin-tip  gap  of  0.35  mm)  shows  the  poorest  per¬ 
formance  among  the  set  of  these  tubes  with  a  fin  density  of 
1.02  fins/mm.  As  visually  observed  during  experimental 
runs,  this  tube  ejected  a  considerable  portion  of  vapor 
bubbles  radially  along  tube  periphery.  This  premature 
departure  of  vapor  bubbles  results  in  a  lower  performance 
by  producing  a  tl  cker  liquid  film  at  the  hot  surface. 

In  comparison  with  the  results  for  tubes  5,  6  and  7,  the 
curves  representing  tubes  2,  3  and  4  show  that  for  a  given 
fin-tip  gap,  the  performance  has  decreased  with  decrease  in 
fin  density.  This  observation  is  believed  to  be  due  to  the 
presence  of  a  greater  number  of  channels  (for  a  given  length 
of  tube)  in  the  higher-fin-density  tubes.  Again,  for  a  fin 
density  of  0.75  fins/mm,  the  0.25  mm  fin-tip  gap  (tube  3) 
resulted  in  the  best  boiling  performance.  The  same 


Figure  26.  Boiling  Performance  of  Gewa-T  Tube 
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Figure  23.  Boiling  Performance  of  Gewa-T  Tube 
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Twa*  T ina/  <8-  late  > 

r  1  d  -  r  <  9  > 

Tv*  T ( 8 ) 

Rho=FNRho<  Cf . Tsat  > 

Mu=FMMu<Cf . Tsat > 

K-FNK <  C  f . Tsat) 

Cp*FnCd  <  Cf . Tsat ) 

3eta=FN3eta<Cf.Tsat) 

N 1  =  rtu/Rho 
Alpha=</<Rho*Oo> 

P  r  =  N  i  /A  lpha 
Amp*V  r /2 .031 
0=Vh»Amp 

Tw  =  Tija-0*L0G<D2/0  1  ) /<  2*8  r*<cu-L  > 

ThetaC  =  Tn'"Tsat 
IF  Thetan<0  THEN 
BEED 

OISP  ,tUALL<TSAT  <HIT  CONTINUE  HHFN  READY)' 

PAUSE 

GOTO  2210 

END  IF 

Hbar  = 1 90 

re=  <  Hbar*p/ <  Kcu.*A  )  )  -5*Lc 

Tanh  =  FNTanh<Fi?) 

Theta* T6etab*7ann/Fe  _  ^ 

<x*<  A  .  81  *6eta*Tneta£>*0o  ?•»  '  anri/ <  Fe»N  1  »A  [pha) ) 
Yy=< i m ,559/Pr >  <9/i6)>  (8/27) 

Hbarc=</Do* < . 6t . 387*yx/Yy >  2 
IF  ABS<  <  Hbar  -  Hbar  c )  /<-  bare  )  >  .  00 1  THEU 
Hbar*  <  Hbar-i-Hbarc )  *  .  b 
GOTO  2720 


.  *  5666  7 
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2'?')  END  'S 

2929  0  I  m  <  Hbar  -P**  cu-A  )  .  5*  r8e  tao-*Tann 

2830  Qc-0-2-01 

j:340  As*PI-02-L 
2850  Qap-Oc/As 

2880  IF  low*  i  THEM  PRINT  USING  " ’  OX  .  DD . 2X  .  ID  .  DD  .  2X  .  2<  l .  aDE  .  2X )  .  0D  .  DD"  :  J  .  ‘m .  Hba 

r . Gao . The tao 

2870  IF  Iow-2  'HEN 

2880  PRINT  USING  "’OX, '"’Data  set  number  *  "".DD" :  : 

2880  PRINT  USING  " 1  OX  .  ""Pos  1 1  ion  numoer  :  1  2  3  4  5  8 

'  3 . 

2900  PRINT  USING  "1  OX  .  "',Temperature  :  ""  .  3 <  DD  .  DD .  i  X  > "  :  !  <  0 )  .  T<  1  )  .  T ( 2  >  .  T <  3  )  .  ' 

<4>.T<5>.'<6).T<7> 

2’R'  0  PRINT  USING  "I  OX .  ""T I  iqd  Ivapr  Tua  Hbar  Qdp  'hetao . 

2920  PRINT  USING  "3X , 3(3D . DD . 2X) ,2<  Z , 4DE .2X ) . 2D . DD" : T Id . Tv . Twa. Hbar . Qap . 'ne tao 

2930  END  IF 

2940  IF  Itn-i  THEN 
2950  SEEP 

2980  INPUT  "OK  T0  STORE  rHIS  DATA  SET  « ! «Y  . 0-N ) 9" . Ok 
2970  END  IF 

2980  IF  3k =1  OR  Ira-2  THEN  J-j+t 

2990  IF  Qk-1  AND  Im»  1  THEN  OUTPUT  4iF  i  ie  :  Emf  <  * )  .  Vr  .  Oh 
3000  IF  Im-2  OR  0k  - 1  THEN  OUTPUT  «P lot :Odp . The tab 
3010  IF  Ira-’  THEN 
3020  3EEP 

3030  INPUT  "HILL  THERE  3E  ANOTHER  RUN  < I = / . 0-N> ?" , Go_:>n 
3040  Nrun-J 

3050  IF  Go_on<>  1  THEN  3H0 
3080  FIAT’  300 

3070  IF  3o_on=  I  'HEN  Repeat 
3080  E  i_  3  E 

3090  IF  jXNran *1  THEN  Repeat 
3100  END  IF 
3  1  1  0  IF  Ira- 1  THEN 
3120  .  BEEP 
3 ' 30  PRINT 

3140  PRINT  USING  " ! OX . ""MOTE :  "".ZZ.""  data  runs  'jere  stored  m  file 


1  .  D_P  i  le'S 
3 '50  END  IF 
7180  8EED 
3170  PRINT 

3180  PRINT  USING  "’OX. ""MOTE:  "".ZZ.""  X-T  pairs  uere 


, 1 0A" :  J- 1  files 
3130  ASSIGN  -File  TO  - 
3200  ASSIGN  •iPiot  TO  - 
32 ’0  SEEP 

3220  INPUT  "LIKE  T0  PLOT  DATA  < I-Y.0-N)?" 

3230  IF  Ok-1  THEN 

3240  CALL  Plot 

3250  END  IF 

3280  SUBEND 

32 70  DEF  FMMutCf.T) 

323U  IF  Cf-1  THEN 


3290' 

RANGE:  28  TO  50 

DEG  C 

3300 

Mu* 8 . 36298 19E- 4- 

-  T*< i . 1094609E-S 

33’  il 

END  Ir 

I:  "20 

IF  Cf -2  'HEN 

3330 

4-24’ .8/<  r-i 33. ' 

!  5  > 

3  ’40 

Mu- 2 ■ 4E -5 *  1 0  A 

7350 

end  :r 

.Ok 


566829E-8) 


stored 


in  plot  data  -lie 
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1  ■>-  V  y  I.  v 


Ls.'L-.'-.'vV 


360 

RE "URN  Mu 

3  370 

FNEND 

3380 

DEF  "NCptCr*.  '> 

3390 

IF  0 f  = 1  THEN 

3400! 

RANGE:  0  '0  30  DE‘ 

3410 

Cp'9 .2507273E- ' +r 

3420 

END  Ic 

3430 

IF  Cf-2  THEN 

3440 

Cp-4.21 1 20858- T  *< . 

3450 

END  IF 

3460 

RETURN  Cp*M00 

3470 

FNEND 

3480 

DEF  FNRho< Cf  .  T) 

3490 

IF  Cf - 1  THEN 

3500' 
3510 
3520 
.3530 
3540 
3550 
3560 
3570 
3580 
3590 
3600 
3610 
3620 
3630 
3640! 
3650 
3660 
3670 
■  3630 
3690 
3700 
3/10 
3720 
3730 
3740 
3750 
3760 
3770 
3780 
3  790 
3800 
3810 
3820 
3330 
3840 
3350 
3860 
3370 
3880 
3890 
3900 
3910 
3920 
39.30 
3940 
3950 
3960 


►<  4 . 4236:  E-5+2 . 7 1  428E.-7  *  r  >  > 


RANGE :  30  TQ  80  DEG  C 

Ro»1 ,6207479E+3-T*(2.2186346'-T*2.3578291E-3> 

END  IF 
IF  Cf-2  THEN 

Ro=999 .52946<-T»<  ,01269-'*<5.482513E-3-T*l ,234'47E-5) > 
END  IF 
RETURN  Ro 
FNEND 

DEF  FNPriCf.T) 

Pt*FNCd  < Cf  , T )*FNMu( Cf . T ) /FNK  <  C  f . T ) 

RETURN  Pr 
FNEND 

DEF  FNK(Cf.T) 

IF  Cf - 1  THEN 
RANGE:  30  TG  90  DEG  C 

K=8 .2095238E-2- T*< 2.221 4286E-4+ T-2 . 3809524E-3 > 

END  IF 

IF  Cf *2  THEM 
X  =  <  T  +273 . 1 5  > / 273 . 1 5 

K  =  -  , 92247  +X*<2 . 8395 -X*< 1  .8007-X*< . 52577- . 07344-X >  > > 
END  IF 
RETURN  K 
FNEND 

DEF  FN T anh  <  X ) 

P=EXP<  X ) 

OEXP<  -  X) 

T  anh  =  <  P-Q ) / ( P+Q ) 

RETURN  Tann 
FNEND 

DEF  FNTvsv  <  V ) 

COM  /Cc/  C <  7  > 

T-i'(0) 

FOR  1=1  tq  7 
T  =  T*C<  D-R  '  I 
NEXT  I 

T  =  t*8 .34643E-  1  -T*(2.45‘j7896E-2-t-i  .  3722S83E-4 > 

RETURN  ' 

FNEND 

DEF  RNBeta< Cf . T> 

Rop  =  FMRho<  C  f  .  T  *  .  1  > 

Rom  =  FNRho<  Cf ,  T-  .  1 ) 

8«?ta=-2/(Rop+Rom)*'RoD-Rom)  /  .2 

RETURN  Beta 

FNEND 

DEF  :rNTaat  <  Cf  . 3  ’ 

:  i  = '  o 
*s» 1 :  n 


•  .y.y 


;V0 

”a=  <  "  .  - > -  .  5 

3  y  8  0 

Pf-  =  rNPsa:  <  C f  .  "a) 

3990 

IF  A6S< <P-4c)/P>> .0001  THEM 

4000 

jr  oc>p  'HEM  'h='a 

4010 

IF  3c<P  rHEN  T 1  =  fa 

4020 

GOTO  397C 

4030 

EMD  IF 

4040 

RETURN  Ta 

4050 

"NEND 

4060 

DEF  i- NPsac  <  C  T  .  Tc  ) 

4061 

IF  Cf-i  THEN 

4070 

r=rc*i .3+32+459.6 

4080 

P- >  0  <  33 . 0655-4330 . 38/T-9 . 2635*LGT<  T  >  +2 . 0539E- 3* 

4081 

END  IF 

4082 

IF  Cf =2  THEN 

4083 

P-Pvs  tu  < T  c ) 

4085 

END  IF 

4090 

RETURN  0 

4  1  00 

FNEND 

41  1  0 

SUB  Plot 

4  ’  20 

DIN  C ( 9 ) 

4130 

INTEGER  Ii 

4  1  40 

PRINTER  IS  1 

4150 

Idv-  1 

4160 

3EEP 

4170 

INPUT  "LIKE  DEFAULT  VALUES  FOR  PLOT  <1-Y. 

0-N)d". 

4180 

BEEP 

4  190 

PRINT  USING  -4X.— Select  Option; . 

4200 

PRINT  USING  "4X.""0  a  versus  aelta-i 

4210 

PRINT  USING  " 4X , 1  "  1  .6  versus  deira-’  ' 

4220 

PRINT  USING  "4X,""2  h  versus  a . 

4230 

PRINT  USING  "4X.""3  h-ratio  versus  delta 

_  J-  «•  •«  •• 

4240 

INPUT  Qpo 

4250 

IF  0po- 3  THEM 

4260 

BEEP 

4270 

INPUT  "SELECT  TUBE  DIAMETER  <0*. 75. 1-1.0 

IN)" . I t 

4280 

END  IF 

4290 

PRINTER  IS  705 

4300 

IF  IdvOl  THEN 

4310 

BEEP 

4320 

INPUT  "ENTER  NUMBER  OF  CYCLES  FOR  X-AXIS" 

.  Ox 

4330 

BEED 

4  340 

INPUT  "ENTER  NUMBER  OF  CTCLES  FOR  Y-AXIS" 

.Cv 

4350 

8EED 

4  360 

INPUT  "ENTER  MIN  X-VPLUE  (MULTIPLE  OF  '0) 

” . Xmi  n 

43  70 

BEEr 

4  380 

INPUT  "ENTER  MIN  Y-VALUE  (MULTIPLE  OF  10) 

** .  Yfn  i  n 

4390 

ELSE 

4  4  Q  0 

0  y  *  3 

441  n 

if  Gdo=o  them 

4420 

Gx  *  2 

44  30 

X»n  t  n  =>  1 

uuai) 

Yni n* 1 000 

4450 

END  I- 

4460 

IF  Ooo-i  THEN 

4470 

C  x  > .; 

4480 

<mn=' 

4  490 

f  n  l  n  *  *  1 .!  f 

•■500 

end  :r 

46  '  II 

IF  0do»J  'HEN 

-520 

'I  X  '  J 
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5 '40  IF  j» :  'HEN  PRINT  "TL  i  S’ 

5 : 50  IF  J>i  THEM  PRINT  "TL  0 

5  i  60  Ya=Yat-»J 
SI  70  Y-LG  I  (  Ya/Tm i n  “Sf  y 
5130  PRINT  "PD  PA  I  00 . " . Y . "YT" 

5190  NEXT  J 

5200  NEXT  I 

52 10  3RINT  "PA  100. ! UO  PET 
5220  PRINT  "PA  0.-2  SR  1.5.2" 

5230  Ii=LGT(Xnun) 

5240  FQR  1-1  TO  Cx+1 
5250  Xa*Xmi n»1 0 ' < I -  1 > 

5260  X*LGT < Xa/Xmi n ) *Sf x 
5270  PRINT  "PA” : X . "  .  0  ;  " 

5290  IF  1 1  >-0  THEN  PRINT  "CP  -2.-2-.LB10  :PR  -2 .2 :  LB"  ;  Ii :  "" 

5290  IF  1 1 <0  THEN  PRINT  "CP  -2 . -2 ; LB  1 0 :PR  0 .2 :LB" :  i i : "" 

5300  Ii-I i+i 
5310  NEXT  ! 

5320  PRINT  "PM  PA  0.0" 

5330  I  l  =LG' <  ym i n ) 

5340  YTO- 10 

5350  FOR  1=1  TO  Cy  + 1 

5360  Ya*ynun»lO  <  I  - 1  > 

5370  Y>LGT<Ya/Yrrun)*Sfy 
5380  PRINT  "PA  0 . "  ;  Y  ,  "" 

5390  PRINT  "CP  -4  .  -  . 25 : LB1 0 : PR  -2 .2 :LB" : 1 1 : "" 

5400  1 1 -I i*l 

5410  NEXT  I 
5420  IF  IdvOt  THEN 
5430  BEEP 

5440  INPUT  "ENTER  X-LABEL" . X  labels 
5450  BEEP 

5*160  INPUT  "ENTER  Y-LA8EL" . Y lape LS 

5470  ELSE 

5480  IF  Opo-O  THEN 

5490  X  labe  IS  =  "Tij-T-3at  OO" 

5500  Y 1 abe 1 3*"q  <U/m  2>" 

5510  END  IF 

5520  IF  Opo-i  THEN 

5530  X labe IS=" Tu- Tgat  <K)" 

5540  YlaDel'S  =  "6  <W/m  2.K)" 

5550  END  IF 

5560  IF  flpo-2  THEN 

5570  Xlabel3=,"a  (W/m  2>” 

5580  Y  laoe  l'5  =  "h  <W/m  2.K>" 

5590  END  IF 

5600  IF  |jpo  =  3  THEN 

5610  X  labe  15  =  "Tij- Tsat  OO" 

5620  YLaDei'5  =  "h(enh>/h<smooth>" 

5630  END  IF 
5540  END  IF 

5650  PRINT  "SR  I  5.2;PU  DA  50. -'6  CP"  : -I  EN< X  labe  l  5> /2 : "0 :LB" : X  labels 
5660  PRINT  "PA  -14.50  CP  0  " : -LEN ( Y labe ! S > 72*5/6 : "0 1  0 . 1 : lB" : Y labe IS 
5670  PRINT  "CP  0.0  DI" 

5680  Reoeat : 

5690  BEEP 

5700  INPUT  "RANT  tq  PLOT  DATA  FROM  a  RILE  (  '  *Y .  0  =  NO"  .Ok 
5710  IF  01- =  i  THEN 
-720  BEEP 

5730  INPUT  "ENTER  7 HE  NAME  OF  ruE  DATA  PILE" . D_‘ i leS 
5740  ASSIGN  *F i  le  T0  D_-L.eS 


5750  BEE? 

5750  BEEP 

5770  INPUT  -ENTER  tHE  BEGINNING  RUN  NUMBER". Md 
5780  BEEP 

5790  INPUT  "ENTER  the  NUMBER  OF  X-Y  PAIRS  STORED” . Noai r 

5800  BEEP 

5310  PRINTER  IS  ' 

5320  PRINT  USING  "AX . ""Select  a  symool : . 

5830  PRINT  USING  "4X.""1  Star  2  Plus  sign . 

5840  PRINT  USING  "4X.""3  Circle  4  Square . 

5850  PRINT  USING  "4X.""5  Rombus . 

5860  PRINT  USING  "4X.""S  R ignt-s ide-up  triangle . 

5870  PRINT  USING  "4X.""7  Up-s  lde-aoun  triangle . 

5880  INPUT  Sym 
5890  PRINTER  IS  705 
5900  PRINT  "PU  DI" 

5910  IF  Sym-1  THEN  PRINT  "SM*" 

5920  IF  Sym-2  THEN  PRINT  “SM  +  " 

5930  IF  Sym  =  3  THEN  PRINT  "SMo" 

5940  IF  Md>  t  THEN 
5950  FOR  1=1  TO  (Md-1) 

5960  ENTER  «File:Ya.Xa 

5970  NEXT  I 

5980  END  IF 

5990  FOR  1=1  TO  Npairs 

6000  ENTER  iFile:Ya.Xa 

6010  IF  0po= 1  THEN  Ya=Ya/Xa 

6020  IF  0po=2  THEN 

6030  Q=Ya 

6040  Ya=Ya/Xa 

6050  Xa=Q 

6060  END  IF 

5070  IF  Gpo=3  THEN  Ya=Ya/FNHsmooth(Xa. ltd) 

6080  X  =  LGT  <Xa/Xmin)»Sf x 
6090  Y=LGT(Ya/Ymin)*Sfy 
6100  IF  Sym>3  THEN  PRINT  "SM" 

6110  IF  Sym<4  THEN  PRINT  "SR  1.4. 2. 4" 

6120  PRINT  "PA"  ,  X  .  Y  , 

6130  IF  Sym>3  THEN  PRINT  "SR  1.2. I. 6" 

6140  IF  Sym  =  4  THFN  PRINT  "IJC2 . 4 . 99 . 0  . -8  . -4 . 0 . 0 . 3 . 4 , 0  :  " 

6150  IF  Sym=5  THEN  PRINT  "UC3 . 0 . 99 . - 3 . -6 . - 3 . 6 . 3 .6 . 3 .  -6 

6160  IF  Sym=6  THEN  PRINT  "IJC0 . 5 . 3 . 99 . 3  . -8  .  -6 . 0 . 3 . 8  : " 

6170  IF  Sym=  7  THEN  PRINT  "UC0 . -5 . 3 . 99 . - 3 . 3 . 6 . 0 . - 3 . -8 : " 
6130  NEXT  I 
6190  BEEP 

6200  ASSIGN  i'File  TO  * 

6210  GOTO  5690 
6220  END  IF 
6230  PRINT  "PU  SM" 

6240  BEEP 

5250  INPUT  "WANT  TO  PLOT  A  POLYNOMIAL  (1  * Y . 0=N ) r' . Go_o 
6260  IF  Go_on= 1  THEN 
6270  BEEP 

6280  INPUT  "ENTER  LOWER  AND  UPPER  X-LIM [ TS" . X  1 1 .  X 1  u 

6290  FOR  Xx  =  0  TQ  Cx  STEP  Cx/200 

6300  Xa=Xmin*IO  Xx 

6310  IF  Xa<X  1  1  OR  Xa>XIu.  THEN  6380 

6320  Ya=FMPo 1 y <  Xa) 

6330  Y=LGT < Ya/ vm i n ) *S f y 
6341)  X=I.GT  (  Xa/'Xm  in )  »Sf  x 
6350  IF  Y < 0  THEM  Y  =  0 
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■jjbO  O-UO  HEN  iilDTC  bo8U 

o370  PRINT  "PA”.x.r ."PD" 

6 ISO  NEXT  < x 
6390  END  Ic 

6400  PRINT  "PU  PA  0.0  SPO" 

SA 1 0  SUBEND 

6420  DEF  FNHsmootni X . I td > 

5430  Hs^ENPo  ly ( X  > /X 

6440  IF  I  td-  1  THEN  Hs=Hs*. 83347 

6450  RETURN  As 

6460  FNEND 

6470  DEF  FNPoly(X) 

6480  Poly=-4. 412371 8E+2-X*<6 .812391 7EJ-2-X«3 . 741 6863E+2) 

6490  RETURN  Poly 
6501  FNEND 

6510  DEF  PNPvst< Tsteam) 

6515  DIM  K<8> 

5520  DATA  -7.691234564.-25.08023696.-168. I  706546 . 64 . 23285504 .- I  13.9646225 
6525  DATA  4 . 1 67 1 1 732 . 20 . 9750676 . 1  . E9 , 6 
6530  READ  < <  *  > 

6535  T- < Ts team+273 . 1 5 ) /647 . 3 

6540  Sum*0 

6545  FOR  N=0  TO  4 

6550  Som^Sum+K ( N)*< I -T )  < N+ I > 

6555  NEXT  " 

6560  8r-Sum/<  T*<  1+K<5>*<1-T>+K<6>*<!-t>  2)>-<1-T)/(K(7)»(t-T>  '2K<8>  > 

6565  Pr -EXP <  Br ) 

6570  P-22120000*Pr 

6575  RETURN  ° 

6580  FNEND 
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APPENDIX  B 


SAMPLE  CALCULATIONS 

Data  set  number  fifteen  of  file  GT6H29  corresponding  to 
tube  number  6  was  chosen  for  the  sample  calculation. 


TEST-SECTION  DIMENSIONS 


O.D. 

HF 

fins/mm  - 

inter-fin  gap  = 
D-, 


U 


25.30  (mm) 
1.05  (mm) 
1.02 

0.25  (mm) 
23.10  (mm) 
17.86  (mm) 
22.22  (mm) 
20.32  (mm) 
31.75  (mm) 
50.80  (mm) 


B.  MEASURED  PARAMETERS 


RES 


=  70.0  (volts) 

=  4.8  (volts ) 

=  2.031  (ohms) 

=  383  (W/m.K) 

=  54.23  (°C) 

=  54.29  (°C) 

=  53.91  (°C) 


91 


T4 

=  53.65  CC) 

1  1  * 

T5 

=  53.40  CC) 

3 

T6 

=  53.31  CC) 

- 

T7 

=  53.22  CC) 

« 

T8 

=  53.54  (°C) 

' 

tsat 

=  47.51  (°C) 

C.  PROPERTIES  OF  FREON-113  AT  SATURATION  TEMPERATURE 

. 

U 

=  4.948x  10  4  (kg/m.s) 

P 

=  1510.0  (kg/m3) 

CP 

=  0.973  (kJ/kg.K) 

k 

=  0.0715  (W/m.K) 

:  .  ■ 

D.  HEAT 

-FLUX  CALCULATION 

The 

heat  transfer  rate  from  the  cartridge  heater: 

mi -■'■M 

qh 

=  VI  (W)  and  I  =  VRES/R  (amps) 

I  = 

4.8/2.031  =  2.363  (amps) 

•-  -  is 

<1 

qh 

=  70.0x2.363  =  165.410  (WO 

The 

cross-sectional  area  of  machined  ends:' 

; 

AC 

=  tt/4 . 0  (D2  -  D2)  (m2) 

O  1 

Ac 

-  y  “"6  2 

=  tt/4 . 0  (0 . 022222  -  0.02032^)  =  63.481  x  10  (m  ) 

92 

The  tube  outside  wall  perimeter: 

P  =  ttD  =  3.1416  x  0.02222  =  6.981  x  lo“2  (m) 

o 

The  corrected  length  of  ends: 

L_  =  Ltt  +  t/2.0  t  =  D  -  D. 

CL  o  1 

Lc  =  0.03175  +  (0.02222  -  0.02032)/2.0  =  32.7  xlO'3  (m) 

Other  thermophysical  properties: 

v  =  u/p  =  4 .948  x  10_4/1510 .0  =  3.277  x  10‘7  (m2/s) 

a  =  k/pCp  =  0.0715/(1510.0x973)  =  4.866  x 10~8  (m2/s) 

1  Ap  1  1509.77  -  1510.26, 

“  c  AT  "  1510  47.61  -  47.41  ' 

=  0.00162  ( 1/K) 

Pr  =  v/u  =  3 . 277  x io_7/4 . 866  x io-8  =  6.734 

Temperature  at  the  base  of  the  T-shaped  fins: 

ln(D  /D  ) 

T  =  T  -  O  - - - — 

B  AVG  “^H  2  7T  L  kc 
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0B  =  TB_TSAT  =  53.342  -  47.51  =  5.832  (K) 


Average  heat- trans fer  coefficient  at  the  unenhanced 
ends  : 


vTn  1/2 


~  '  0.60  +  0.  387 


gS^eBTanh{;Ep-)  ■ 

T  .  hP  . 1/2 

Lc  -* 

[1  +  (0 . 559/Pr) 9/16] 8/27 


h  =  187.49  (W/m  .K) 


Heat- transfer  rate  through  unenhanced  ends 


vh  P  kc  Ac  0g  Tanh  mLc 


(  hP  } 1/2  _  (187.489  x 6.981  x IQ  2  ^ 1/2 

kCAC  0.0715  x  63.481  x  10_6 


=  1.698  x 10 


Qf  =  3.29  (W) 


QLOSS  =  2qf  =  6*579  (w) 


Heat  flux  through  Gewa-T  surface: 


Q  =  165.41  -  6.579  =  158.83  (W) 

q  =  Q/Ag 

Ag  =  ~r  * 23  .  1  >  lO"3  x  50 . 8  x  lo“ 3  =  3 . 687  .«  10_3  (m2) 

ICO  Ol  O  9 

q  =  - _  =  43.978  x  10J  (W/m  ) 

3.687  >  10 

The  following  are  the  results  obtained  from  the  computer 
y  running  the  data  reduction  program  (DRP)  (see  Appendix  A) 

q  =  44.466  10  (W/m2) 

0B  =  5.76  (K) 

h  =  187.39  (W/mZ . K) 


The  small  difference  in  results  obey  the  inaccuracy  of 
and  calculations  with  the  use  of  about  three  significant 
igits  only. 
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END 


2/2 

NL 


m 

Ui 


HIM 

12.2 

2.0 
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APPENDIX  C 


UNCERTAINTY  ANALYSIS 

The  same  data  set  (set  No.  15  of  file  GT6H29)  that  was 
used  for  the  sample  calculation  was  chosen  for  the  uncertainty 
analysis;  therefore  dimensions  of  test  section,  measured  and 
calculated  parameters  found  in  sample  calculation  were  used 
in  this  analysis.  All  uncertainties  are  presented  as  a 
percentage  of  the  calculated  parameter. 

A.  UNCERTAINTY  IN  SOURCE  HEAT- TRANSFER  RATE 


Q  =  VI  (W) ,  I  =  vres/R  (amp) 

6V  =  6VR£S  =  0.05  (volts) 

where : 


5  = 

51  = 

5Qh 

QH  = 
6QH 

qh 


uncertainty  in  measurement  and  calculation 


<5VRES/R  =  °-05/2-031  =  0.025  (amp) 


III  2  r  T  2  1/2 

nil)  <il)  i 

V  I 


,0.05,2  ,  ,0.025,2, 

~~To~  +  (27363)  J 


1/2 


1.06  percent 


UNCERTAINTY  IN  SURFACE  AREA 


Ag  =  tt  D2  L  6D2  =  6l  =  0.1  (mm) 


SAg 


[( 


0.1  . 

23. 1J 


+ 


0.1 

50.8 


2  1/2 
)  1 


0.48  percent 


UNCERTAINTY  IN  AT 


AT 


T 

W 


'SAT 


5TW  =  5TSAT  -  0.15  (°C) 


6AT 

AT 


6Tw  2 

[(— -)  + 

AT 


(- 


6T 


SAT, 


2  1/2 

1 


AT 


6  AT 
AT 


ug4l) 2  +  (-  §4l)2i 


1/2 


6.18 


6.18' 


3.43  percent 


UNCERTAINTY  IN  HEAT  FLUX 


q 


H 


~  2qf 

*13 


= 

q 


[( 


6Q. 


qh'2qf 


+  ( 


2<5Q 


Qh-2qf 


F  )2  +  (- 


6A  ?  1/2 

— )  ] 

A  J 

ab 


In  this  case,  Qrr  =  165.41  (W)  and  Q_  =  3.29  (W) 


Assuming  the  same 


Therefore,  Q  =  2.0  percent  of  Q„. 

r  n 

proportion  in  the  uncertainty  for  Q  : 

F 

6Qf  =  0.02  6Qh  =  0.02x1.753  =  0.035  (W) 


^a  = 
q 


r  ,  1.753,  2  ,  ,  2  x  0.035, 2 

158.83  158.83  1 


+ 


(-0.0048) 


1/2 


Aa  = 
q 


1.26  percent 


3.65  percent 
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